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Abstract 

By employing a Boltzmann transport equation and using an energy and size dependent relaxation 
time (r) approximation (RTA), we evaluate self-consistently the thermoelectric figure-of- merit ZT 
of a quantum wire with rectangular cross-section. The inferred ZT shows abrupt enhancement 
in comparison to its counterparts in bulk systems. Still, the estimated ZT for the representative 
Bi2Te3 nanowires and its dependence on wire parameters deviate considerably from those predicted 
by the existing RTA models with a constant r. In addition, we address contribution of the higher 
energy subbands to the transport phenomena, the effect of chemical potential tuning on ZT, and 
correlation of ZT with quantum size effects (QSEs). The obtained results are of general validity 
for a wide class of systems and may prove useful in the ongoing development of the modern 
thermoelectric applications. 

PACS numbers: 73.50.Lw, 73.23.-b, 73.63.Nm 
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I. INTRODUCTION 



Thermoelectricity deals with conversion of heat to electricity and vice versa. Last few 
years, mainly due to the recent advances in material science and fabrication techniques on 
the one hand and a quest for alternative sources of energy generation on the other hand, 
have witnessed a rejuvenation of interest in the thermoelectric phenomenar-^. Performance 
of a thermoelectric device is defined through its dimensionless figure-of-merit ZT = S^aT / n. 
Here, T is the absolute temperature, S stands for the thermopower (the Seebeck coefficient), 
and a and k represent the electrical and thermal conductivities, respectively. At first glance, 
a large ZT might be achievable by maximizing a and S*, and by minimizing k. The problem, 
however, lies in the fact that these quantities are interdependent and each of them is rather 
sensitive to the material properties, to the temperature, and to the size and geometry of the 
underlying building blocks in a miniaturized thermoelectric module. Consequently, often 
one of the two not-completely-independent strategies is exploited: increasing of the power 
factor P = S^a which is largely an electronic property, or decreasing of k, which consists of 
electronic and lattice Kph. contributions by manipulation of the latter. A large P can be 
achieved by tuning of the chemical potential level through doping or external electric fields 
(gating) and by having control over the scattering mechanisms of the charge carriers in the 
device. Reduction in Kph^, instead, can be obtained by exploitation of material systems built 
up from heavy elements and by intentional enhancement of the phonon scattering events 
through introduction of lattice imperfections (impurity atoms, superlattice structures, etc.) 
or by intensifying of the surface scattering of the heat carriers from the specimen boundaries 
in the systems with reduced dimensionality. 

In what follows, assuming a constant Kph., we focus on the charge carrier contributions to 
the thermoelectric properties of a one-dimensional (ID) wire with discrete energy levels. To 
demonstrate the applicability of the present approach, we evaluate ZT in bismuth telluride 
nanowires and discuss the obtained results with an eye on the predictions of the earlier 
models. 
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II. THEORETICAL MODEL 



Current generic models (as compared to the ab initio calculations) for evaluation of 
the thermoelectric figure-of-merit are mostly based on the solution of Boltzmann transport 
equation (BTE) with various degrees of sophistications. An efficient approximation for 
solving BTE can be attained by its linearization and through introduction of a relaxation 
time r, the time period within which the system gains its equilibrium after removal of the 
external stimulus^^. Consequently, the components of the transport tensors of the system 
can be obtained through the relation 

[C^%s = ^ / {-dEfo)r{k){E{k) - C)Vv,dk, (1) 

where a = iZ'^^^ designates the electrical conductivity, S = f^fjny is the Seebeck coefficient, 
and Ke = T^f-C^^^ — ^^-^0—] is the electronic contribution to the thermal conductivity^. 
Above, q is the electrical carrier charge, -E'(k) is the energy- wavevector dispersion relation, 
/iv(k) = Vk-E'(k) defines the velocity operator [h is the reduced Planck constant), and 
u is the valley degeneracy. The equilibrium Fermi-Dirac distribution function is given by 
fo{E) = [l + exp f3{E — , where (3 = l/ksT {ks is the Boltzmann constant) and ( denotes 
the chemical potential. 

The precise form of the dispersion relation above can only be derived through detailed 
band structure calculations and depends strongly on the material properties and on the 
boundary conditions in the system of interest. Similarly, the relaxation time for any specific 
process depends on the density of energy states (DOS) and on the types of the scattering 
mechanisms taken into account. In most calculations, however, one expresses DOS via the 
single-particle eigenenergies and assumes a power-law dependence of the relaxation time on 
energy, t{E) oc E". It was only under such conditions that A. F. loffe arrived at his com- 
pact expressions in terms of the Fermi-Dirac integrals for the evaluation of the transport 
coefficients in Eq. ([1]) above^^. Despite this, such basic assumptions have been widely 
overlooked in the latter evaluations of the ZT values in low-dimensional systems, leading 
ultimately to the estimations of unrealistically high, monotonously size-dependent, ZT's 
in atomically thin nanowires^"— . It is worth reminding that in the limit of ultra narrow 
nanowires, the whole concept of a size-independent DOS and any argument subsequent to 
it will collapse. In such cases, one has to resort to the other proper techniques like, e.g., 
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tight-binding models and check for the vahdity of the obtained results self-consistently. Fur- 
thermore, as shown elsewhere^i, even minor variations in the form of DOS [and subsequent 
changes in t{E)] will affect ZT values dramatically. It is also elucidating to remind that, 
based on a purely mathematical analysis of the functional in Eq. ([T]), Mahan and Sofo 
estimate in2i a universal (that is, independent of size, temperature, and specific material 
properties) finite upper limit for ZT in semiconductor thermoelectric materials. Based on 
these facts and to circumvent the above-mentioned crucial shortcomings in the existing theo- 
retical models for low- dimensional thermoelectric systems, starting with the basic principles, 
below we obtain the transport coefficients and ZT values specifically for a ID nanowire with 
quantized energy levels. The key issue in such a treatment is the derivation of a proper, i. e. 
size and energy dependent, expression for the relaxation time. 

Let us assume that scatterers each with a scattering strength Vq are randomly dis- 
tributed at positions Rj along the wire, V{r) = J2f=i^o^i'^~^j)^ make use of the 
Fermi's golden rule T~^jr = (2n/h)\{i \ V{r) \ f)f to derive the corresponding scattering 
rates. By averaging over the configuration of the scattering centers^^, we arrive at 

where A^J^, = (2 + Smm')^^ + Snn') and g = N/Q. Here, Q=wtL is the volume of the wire; 
w and t are the lateral dimensions and L stands for the wire length. The energy density of 
states per unit volume is given by 

1 [7^^QiE-E,^ 



Here, is the component of the charge carrier mass along the wire axis, s is the spin 
degeneracy, and G designates the Heaviside step function. Above, use has been made of 
the fact that the single-particle wave functions and the corresponding eigenenergies of a 
carrier confined in lateral dimensions x and y and traveling freely along the z-axis are given 
by '^ij{k) = (2/v^) sm(iTTx/w)sm(jTry/t)exp(ikz) and E=Eij + (h'^k'^ /2mz), respectively; 
Eij = [fx^Ti"^ /2){j? /w'^nix+j'^ /t^my) are the subband energies, i, j G N, and k is the wavenum- 
ber. 

In Eq. ([2]) above, as in most cases the density and strength of scatterers are unknown 
and may vary from one individual wire to other, one can make use of the definition of carrier 
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mobility ^=q (r) /m* and account for the product qVq througb^^ 

1 sv VL fim* 



4 wt qlT? V 2kBT ^ (r„„) 



(4) 



mn^m' n' 



where the hfetime of the state {m,n) is defined through the relation = ^rn'n''''nm- 
and its expectation value as (rmn) = fTmn{£)£g{£)dsfde/ Jeg{e)defde. Here, e = f3E is 
the reduced energy and m* is the effective mass (for a ID wire aligned along the 2;-axis, 
m* = rriz). Now, substituting for the corresponding expressions and making a coordinate 
transformation where e and the reduced chemical potential (* = (3( (for electrons) are 
measured from the bottom of the conduction band^, the transport coefficients in Eq. ([T]) 
can be expressed as 



(;) su Awn m* Irriz q 
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where J"- - = TUC) and 
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(6) 



Total contribution from all the subbands will now be = Ylij^ij • (The contribution of 
the holes can be accounted for analogously.) 

This equation has to be compared to the corresponding ones based on a size-independent 
power-low RTA approach introduced originally in^ and exploited later by Dresselhaus and 
coworkers in their pivotal^ and subsequent studies^i^ on low- dimensional thermoelectric 
systems. 



III. RESULTS AND THEIR DISCUSSION 



Next, we evaluate ZT for a bismuth telluride nanowire at T = 70 K and at 300 K. To 
make the comparison between the predictions of the present model and those of the previous 
studies^"— more transparent, we use the same material parameters, = 0.32, rriy = 0.08, 
and m* = rriz = 0.02 for the electronic effective masses (all in units of the free electron mass), 
and take the same value for the mobility of electrons, fiz = 1200 cm^V~^s~^. The lattice 
thermal conductivity of bulk Bi2Te3 is k^JJ"^ ~ 1.5 W/Km. The phonon confinement effects 
become considerable only if the lateral dimensions of the wire, w and t, are comparable in 
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size to the phonon mean free path Xph.. As for bismuth telluride nanowires Xph. ~ 1 nn>i^, 
for wires with w,t ^ Xph,, one can safely assume Kph. ~ k^^—- Figure 1 shows dependence 
of the electronic contribution to the figure-of-merit of a wire with square cross-section on its 
size w and on Q* . Here, in contrast to the constant-r model of Hicks and Dresselhaus (HD) 
(see, the Appendix), dependence of ZT on wire size is non- monotonous and after reaching 
its optimum at a certain width and chemical potential, falls starkly to zero. This is a direct 
consequence of disappearance of electronic density of states for ultra thin wires. The most 
distinct differences between the two models are as follow: the peak values of the a, Kg, and 
ZT in the present model show a shift toward larger wire thicknesses, ~ 10 nm, as compared 
to the plain w~^-dependence of the corresponding quantities predicted by the HD model; the 
thermopower S given here is temperature dependent and reaches the overall high values of 
the HD model only in very narrow wires; in contrast to the HD model, here the position of 
optimum ZT in (^w-plane has a visible dependence on temperature and, as temperature falls 
down, shifts toward higher thicknesses; the asymptotic values of the transport coefficients 
for a model square wire with i/; = 50 nm given by the HD model are unrealistically smaller 
than the those obtained here (Table 1; see also the Table in Appendix). 

To illustrate the effect of subband contributions and that of QSEs, the calculations were 
repeated for the first three subbands, up to (i,j) = (1,3) (see. Fig. 2). Table 1 summarizes 
the results. Here, especially at T = 70 K, one can distinguish two clear maxima for a, 
and S which arise from QSEs^^i^^ and which are smeared out in ZT. This observation is 
rather general for other subbands too and is a consequence of the fact that a and have 
similar functional dependencies on the wire size and the QSEs distinguishable independently 
in each of them (and to a lesser extent in S) are compensated in ZT by one another. Another 
interesting observation is the abrupt enhancement of Kg at thinner wire diameters and in the 
region where C,* approaches its optimum. Here, noticeably, k^. ^ np^ . Alike, the HD model 
predicts an even stronger enhancement of in thinner wires (Appendix). Dependence of 
the thermal conductivity on dopant concentration (cf. tuning of C), leading ultimately to 
the takeover of Uph. by the electronic contribution Kg, was observed in experiments with 
bismuth-antimony alloys some years ago2^. 
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FIG. 1: (a) Dependence of the transport properties and the thermoelectric figure-of-merit ZT 
(only the electronic contribution) on wire width w and on the chemical potential (" (measured in 
units of thermal energy k^T^ for a square Bi2Te3 nanowire at T = 300 K. The valley degeneracy is 
1/ = 6. Only the first subband contribution is taken into account, (b) T = 70 K. Notice the many 
fold decreases of the peak values and their shift toward larger wire thicknesses. 
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FIG. 2: (a) Dependence of the transport coefficients on the wire size and on the chemical potential 
at T = 300 K. Same parameters as in Fig. 1, except that here {m,n) = (1,3). The clearly visible 
bumps in a and Kg are manifestations of the quantum size effects, (b) T = 70 K. 
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TABLE I: Values of the transport coefficients for a square 50-nm-tliick Bi2Te3 wire (see Text) 
with valley degeneracy i/ = 6 at a fixed reduced chemical potential = —3. The indices 
correspond to the subband numbers up to which the electronic (vs. holes) contributions are taken 
into account (Kp/^. = 1.5 W/Km). HD stands for the one-band constant-r model employed in 
previous studie o^^'^^" — . For more details, see the Table in Appendix. 

IV. SUMMARY 

The approach described here, which is based on the formulation of the scattering rates 
of the charge carriers in terms of an energy and size dependent expression for the density of 
states, can be readily extended to address transport and thermoelectric properties of a wide 
range of quantum systems subject to more sophisticated forms of the relevant DOS's. The 
model also holds potential to account for the non-diagonal contributions of the transport 
tensor elements to the kinetic properties. It can equally be exploited in the metallic or 
semiconducting regimes and it considers inherently the important issue of quantum size 
effects in low-dimensional structures. 

This work is dedicated to the memory of Amirkhan Qezelli, the uncle, the childhood 
friend. 
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V. APPENDIX 



The complementary Table 2 outlines contributions of the higher energy subbands to 
the transport properties and to ZT in a 50-nm-thick Bi2Te3 nanowire with effective mass 
components m^. = 0.32, vriy = 0.08, and with = 0.02 along the wire axis (all in units of 
the free electron mass). We also summarize the main results of the constant-r model and 
demonstrate its predictions in the appended Figure 3. 



Constant relaxation-time model 

Derivation of the transport properties in bulk systems, based on a constant relaxation- 
time approximation, was originally given in^^ and later discussed in detail by Nolas et al. 
ini^. Later, extension of these results to the case of one-dimensional nanowires was obtained 
by Dresselhaus and coworkers^^ii^. A more recent account is given iri3. Below, we summarize 
the results. 

The electronic (vs. holes) contribution to the transport coefficients cr = (electrical 
conductivity), S = -(l/gT)£(i)/£(o) (the Seebeck coefficient), and = {l/q^T)[C^^^ - 
(C^^^y / C^^^] (the electronic thermal conductivity) are given by the following expressions. 
Here, q is the magnitude of the elementary charge. 



where 



and 



= v/m*/ie (7) 



F.^F.(n^l 5;^. (8) 







Above, e = E/ksT is the reduced energy, m* and fig stand for the electron effective mass 
and mobility along the wire axis, u represents valley degeneracy, and (* = C^/ksT is the 
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FIG. 3: (a) Thermoelectric figure-of-merit ZT for a square Bi2Te3 nanowire at T = 300 K obtained 
by the constant-r model. The chemical potential C is measured in units of thermal energy k^T . 
The effective mass values are the same as in the Text {rrix = 0.32, my = 0.08, and niz = 0.02 along 
the wire axis), (b) T = 70 K. Notice that the predicted enhancement of the transport properties 
here emerge only in ultra thin wires {w ^ 1 nm). This has to be compared to the predictions of 
our model presented in the Text. 
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reduced chemical potential of electrons measured from the bottom of the conduction band. 
The above expressions can be further simplified to: 

' '1/2 



q \R 



1 2u f2kBTy\,^^ 5 9F,% 



The dimensionless figure-of- merit, ZT = o'S'^T/{ne + i^ph.), can now be written as: 



1 f 3Fi/2 p 
B + 2^3/2 



B ' 2F_i/2 

with 

TTW'^ \ J qKph. 

The most distinct difference of the above expressions with the those presented in the Text is 
lack of oscillatory dependence of the transport coefficients on the wire size and on the quan- 
tization of the energy levels, i.e., absence of the quantum size effects. Also, it is noticeable 
that here the Seebeck coefficient S is wholly independent of size and the temperature, and 
a and have only a bare monotonic dependence of the form oc on the wire width. 

Such dependencies lead ultimately to unrealistically large values for a, and, eventually, 
ZT in ultra thin wires. Figure [3] presents the obtained results for a Bi2Te3 nanowire with 
the same parameters as given in the Text. 

This work is dedicated to the memory of Amirkhan Qezelli, the uncle, the childhood 
friend. 
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TABLE II: Electrical conductivity a, electronic thermal conductivity Kg, the Seebeck coefficient 
S, and the thermoelectric figure-of-merit ZT for a quadratic Bi2Te3 nanowire (aligned along the 
z-axis) with w = 50 nm at a fixed chemical potential = —3 (measured from the bottom of 
the conduction band). The indices (m,n) correspond to the subband numbers up to which the 
electronic (vs. holes) contributions are taken into account (/tp/j =1.5 W/Km). HD represents the 
one-band constant-r model^^i"— . Note that with valley degeneracy = 6, a and Kg become simply 
six times larger while S remains intact. Dependence of ZT on u is of the form ZT oc (Kg -|- Hph./^)~^ 
and its values for v = 1 and v = Q are given explicitly. 
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